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Abstract 

The diffusion of glycerol molecules decreases with decreasing temperature as its viscosity in- 
creases in a manner simply described by the Stokes-Einstein(SE) relation. Approaching the glass 
transition, this relation breaks down as it does with a number of other pure liquid glass formers. 
We have measured the diffusion coefficient for binary mixtures of glycerol and water and find that 
the Stokes-Einstein relation is restored with increasing water concentration. Our comparison with 
theory suggests that addition of water postpones the formation of frustration domains pj]. 
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Most liquids become solids after they are cooled through their melting temperature, T^, 
following a first order thermodynamic transition. A notable exception is a class of molecular 
fluids known as glass formers which can stay in a metastable state as supercooled liquids to 
temperatures substantially below their equilibrium melting transition. The thermodynamic 
properties of the supercooled liquid are a continuous extrapolation of the liquid state from 
above T^. On further cooling the dynamical and transport behavior, viscosity, diffusion, 
and motional relaxation times, exhibit such dramatic slow-down over a narrow temperature 
range, that there appears to be a transition to a rigid but amorphous solid state at a 
characteristic temperature, Tg, called the glass transition|2]. 

According to the Stokes-Einstein (SE) relation the product of viscosity and diffusion 
coefficients is simply proportional to the absolute temperature j"" 



D = ^, (1) 

In Eq.l, D is the translational diffusion coefficient, rj is the viscosity, ks, the Boltzmann 
constant, and the hydrodynamic radius. The basis for the SE-relation is that the elemen- 
tary units of the liquid are just Brownian particles and that the system is homogeneous. 
It is remarkable that many glass formers still follow this behavior even as the relaxation 
processes slow-down. Ultimately, at low enough temperature, it fails, at least for a number 



of the more fragile 



J] glass forming liquids. 



Kivelson[l| et al. have developed a theory for frustation limited domains (FLD) that 
form at a temperature T* above the melting temperature T^. Although these domains 
are energetically favorable on a local scale their existence affects diffusional and rotational 
relaxation differently giving rise to breakdown of the SE-relation. In this letter we present 
diffusion measurements of four aqueous solutions of glycerol and compare with viscosity 
measurements to show that the SE-relation fails for glycerol; but that this relation can be 
restored at low temperatures through addition of water. With the FLD theory we extract 
the temperatures for domain formation, T*, for the four solutions and we show they are 
consistent with the temperatures where we first observe deviation from the SE-relation, 
providing a direct confirmation of these theoretical ideas. We have chosen binary solutions 
to study glass formation since they allow systematic and continuous control of the essen- 
tial molecular interactions involved, an approach that may also prove useful for numerical 
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simulations. 

Among glass formers, glycerol is well studied having an intermediate fragilityj^ m = 
53, glass transition temperature T„ = 190 K and melting temperature = 290 K. It is 

n Q n p n 

strongly resistant to crystallization [g, |7j as are its aqueous solutions [8, |9|, [lOj- For glycerol 
concentrations above 55% they are good single-phase glass formers The concentration 
dependent glass transition temperature has been measured [12I, Q and can be described 
quite well by an empirical formula [l^: Tg = 1901^1 + 136^2 + 27wiW2, where wi and W2 
denote the weight concentration of glycerol and water respectively, such that Tg smoothly 
decreases with increasing water content as is shown in Fig. 2. 

Samples were prepared with 100, 95, 90 and 85% glycerol by weight, using pure glycerol 
(99.9%, Fisher Scientific) and D2O (99.9%, Cambridge Isotopes). The NMR signal, which 
we used for measurement of diffusion, was therefore selectively sensitive to the diffusion of the 
glycerol molecule. The mixtures were loaded into pyrex glass ampoules in the dry nitrogen 
atmosphere of a glove bag. The ampoules were fiame-sealed while the liquid mixture inside 
was frozen with liquid nitrogen. The diffusion measurements were performed by proton 
NMR stimulated-echo sequences in the fringe magnetic field, H = 4.28 T, frequency of 
182 MHz, of a superconductive magnet with a magnetic field gradient G = 42.3 T/m. In 
order to improve signal-to-noise, and hence resolution in diffusion, we used the frequency 
jumping method developed by Sigmund and Halperin[l5|,Q, allowing the near-simultaneous 
collection of data from 17 independent NMR-sample slices, each 110 microns thick. 

The slow-down in dynamics of a glass former is marked by deviation of the temperature 
dependence of diffusion from the Arrhenius behavior that is typical at high temperatures: 
oc e~^/'^^^. The modified temperat 
Tamman-Fulcher (VTF) equation 



rure dependence can usually be represented by the Vogel- 



13,1181, [l9j which, for the diffusion coefficient D, is given 



by. 



D = A^Vfexpi--^). (2) 



Here Aq, B, and Tq are phenomenological constants determined from experiment. At still 
lower temperatures, close to Tg, the non-exponential features of relaxation, diffusion, and 
viscosity are often best captured by stretched exponential behavior. A theoretical basis that 
describes the continuous evolution of glass forming dynamics with decreasing temperature 
has not yet been established. Nonetheless, it is clear that the molecular interactions must 
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play a crucial role. 
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FIG. 1: The logarithm of the diffusion coefficient D of glycerol as a function of inverse temperature 
for various aqueous solutions of glycerol. The curves are fits to the VTF-formula. 

The logarithm of the diffusion coefficient for glycerol is plotted as a function of the 
inverse temperature in Fig.l to emphasize deviation from Arrhenius behavior which would 
be a straight line in this figure. The diffusion coefficients shift to the right with increasing 
water content while retaining a temperature dependence similar to that for pure glycerol, at 
least in the high temperature region. But at lower temperatures the diffusion for solutions 
is softened compared to that of pure glycerol. To get a quantitative assessment, we fit the 
experimental diffusion coefficients to the VTF-relation and present the important parameters 



for these fits in Fig. 2 along with T*, determined from the viscosity data[20|, |2l| analyzed 
using the frustration limited domain theory (FLD) and r^, the hydrodynamic radius 
obtained from the SE- relation, Eq.l. 

The VTF-relation, Eq.2, can be derived from a free volume theory {3] • The prefactor Aq 
in this relation is a constant, at most a weak function of water concentration; the activation 
energy, B, was found to be a weak, but systematically decreasing, function of water concen- 
tration; and the temperature, Tq, corresponds to full arrest of the molecular dynamics at a 
temperature where there would be zero free volume. This temperature should be less than, 
but closely related to, the glass transition temperature j^]. A direct comparison between the 
two is given in Fig. 2, showing that Tq scales with Tg. 

The VTF-fits shown in Fig. 1 were performed for diffusion coefficients in the high tempera- 
ture region. They deviate from the experimental values at low temperature where D ^ 10~® 
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FIG. 2: VTF-parameters, the domain formation temperatm'e, T*, and the hydrodynamic radius, 
rg, as a function of glycerol concentration in aqueous solutions. The solid curve is the glass 
transition temperature Tg, taken from the empirical formula given in the text. It closely parallels 
the VTF-parameter Tq. The dashed lines are guides to the eye. 

cm^s~^. In order to understand the failure of the VTF-formuIa at low temperature, we 
return to the free volume theory from which the VTF-relation is derived. In this theory, the 



diffusion of the molecule is associated with the potential cage formed by its neighbors [22], 
such that the molecular interactions are described in a mean-field way. In the case of glyc- 
erol and its aqueous solutions, the dominant molecular interaction is hydrogen bonding 
(H-bonding) and the potential cage should arise from these bonds. Numerous theoretical 
investigations on H-bondin g ha ve been conducted for glycerol, water and other H-bonding 



dominated alcohols 



23 



24 



25 



26| and it has been shown that the detail of H-bonding 



and its effects are complex: one molecule can have different H-bonding states depending 
on the number of bonds; and H-bonding can be distinguished as either intermolecular or 
intramolecular 231], affecting diffusion differently. The life time of the H-bond is also tem- 
perature dependent 23|, l2J] . The more intermolecular H-bonding a molecule has, the more 
of an obstacle there is to diffusion. This intuitive idea is confirmed in calculations by Marti 
et all 



for water and by 
calculations by Chelli et al. 



Matsumoto et al. 



261 ] for methanol. The molecular dynamics 



231 ] and Padro et a/. [24^] indicate that the average number of 
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H-bonds per molecu 
about 6 times longer 



e is higher for glycerol than water, and their life time in glycerol is 



24 [. The average life time of a H-bond in glycerol can be described by 



2J. This 



an Arrhenius equation[23|] but the state of H-bonding also varies with temperature 
complicated situation goes beyond the simple free volume picture [22] and leads to failure of 
the VTF-relation at low temperature where these additional temperature dependences come 
into play. 

Fits to the SE-relation in the high temperature limit give the results for in Fig. 2 with 
a satisfying consistency for different water mixtures but systematically increasing by 25% 
over the range of water concentrations. In terms of H-bonding we argue that a dynamical 
complex of glycerol and water in binary mixtures is the relevant elementary unit, akin to the 
ideal Brownian particle, having a larger effective radius than pure glycerol which increases 
systematically with water content. The increase in the diffusion coefficient of glycerol with 
dilution in aqueous mixtures can be ascribed to the decrease in the effective interaction 
between these Brownian clusters that we associate with the weaker intermolecular H-bonding 
in water. 

In Fig. 3 we show the temperature dependence relative to the glass transition of Drj/T, 
normalized at high t emp erature. The viscosity data for pure glycerol are taken from smooth 
fits of two data sets 20|, l2l|] which are consistent with each other to ~ 3%. The viscosity 



of the aqueous solutions are taken from smooth fits to the data by Segur et. a/. [20]. At 
temperatures T/Tg < 1.62 we resolve a breakdown in the SE-relation for pure glycerol 
beyond the estimated measurement error in diffusion which is represented by the scatter in 
our data. The temperature at which this breakdown is observed moves to lower temperature 
with increasing amounts of water. In fact, for the 85% sample the SE-relation holds over the 
entire range for which there are viscosity data available. The behavior of pure glycerol we 
report here matches qualitatively that of other glass formers where larger deviations from 
the S&relation are observed at temperatures closer to their glass transition. Collectively, 
these results on glass formers are consistent with theoretical models that have recently been 
proposed[27] suggesting failure in the SE-relation irrespective of fragility. 

For pure liquids, such as glycerol, Arrhenius behavior is observed at high temperatures in 
both diffusion and viscosity. So it is natural to expect the SE-relation to hold there. But at 
lower temperature, where VTF-behavior is observed, it is not obvious that the SE-relation 
should remain valid; or if not, what might be the mechanism responsible for its failure. One 
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FIG. 3: Breakdown of the Stokes- Einstein relation through comparison of diffusion and a smoothed 
interpolation to viscosity data 



ycerol. For comparison the diffu- 



se] with the same normalization. 



21[ for aqueous solutions of gl 
sion data of Hrovat and Wade are also shown for pure glycerol 
The validity of the Stokes-Einstein relation is limited to high temperatures. The lower limit of 
temperature in this figure is determined by availability of viscosity data in aqueous mixtures. The 
four arrows on the top axis indicate T*, the domain formation temperature. From left to right 
they are for 85, 90, 95 and 100% glycerol solutions, and are qualitatively consistent with onset for 
breakdown of the SE-relation. 

important aspect of this problem is that the viscous relaxation process is microscopic while 
the measurement of diffusion is macroscopic, meaning that the former is on the molecular 
length scale and the latter on the length scale of tens of microns. So a possible explanation 
for breakdown of the SE-relation could be ascribed to local domain or cluster formation 

29,|30|. 
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rustration limited domains below a 

2a 



2l| using the FLD theory [l| to 



which can lead to dynamical heterogeneity jl|, 

Kivelson et al.\l^ suggested that the breakdown of the SE-relation is due to dynamical 
heterogeneity, resulting from formation and growth of 1 
temperature T*. We have analyzed the viscosity data 
find T* plotted in Fig. 2, and given in Fig. 3 by vertical arrows. The temperature where 
deviation from the SE-relation appears in Fig. 3, is qualitatively consistent with the domain 
formation temperature T*. In particular they both decrease systematically with increasing 
water content, suggesting that the breakdown of the SE-relation is controlled by the birth 
of these frustrated domains. 
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It is now well-established experimentally that the SE-relation fails for certain fragile glass 
formers in their pure form such as ortho-terphenyl l| and tris-napthylbenzene j32 j. and now 
for glycerol, a significantly less fragile material. However, Chang and Sillescu[33|] did not 
find such breakdown behavior for glycerol. This is in contrast to our preliminary workjl6| 
and that which we report here. There are significant differences (~ 30%) in the absolute 
values of the diffusion coefficients for pure glycerol measured by Chang and SillescujsS] as 



compared with other workers 
and Wade 



34 



35l | including the more precise measurements of Hrovat 



341 ] which we have confirmed and extended to lower temperatures; see Fig.3. This 



discrepancy is outside the possible error in our measurements of diffusion and it is most 
likely that there is systematic error in the Chang and Sillescu data that accounts for their 
claim that glycerol obeys the SE-relation over this temperature range. We also note that 
the work of these authors is inconsistent with the FLD theory. 

In conclusion, we have measured the diffusion coefficient of glycerol in aqueous solution. 
We have found that VTF-behavior fails as a phenomenological form for the temperature 
dependence of the diffusion coefficient, but with a range of validity that extends to lower 
temperature with increasing water content. Similarily, we find that the SE-relation breaks 
down for these aqueous mixtures but at a temperature that decreases progressively with in- 
creasing water content. These observations are consistent with expectations from frustration 

n □ 

limited domain theory |27|]. Finally, we remark that the mechanism responsible is most 
likely the heteromolecular hydrogen bonding which postpones the formation of domains and 
inhibits the breakdown of the Stokes-Einstein relation as temperature is decreased. We find 
that the effective intermolecular interaction from hydrogen bonding can be conveniently 
controlled using binary solutions. 

We thank M.D. Ediger for helpful suggestions and acknowledge support from the Mate- 
rials Research Center at Nortwestern University provided by the NSF through grant DMR- 
0076097 and support from the Department of Energy, contract DE-FG02-05ER46248. 
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